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ABSTRACT 



The present invention offers methods and apparatus for 
sonic borehole logging of formation shear using discrete 
frequency measurements. A multipole source trans- 
ducer emits sonic energy at a location in a borehole of 
at least one discrete frequency which induces propaga- 
tion of multipole sonic waves in the borehole; sonic 
energy of the multipole sonic waves is detected at multi- 
ple receiver locations in the borehole by multipole re- 
ceivers to produce a received sonic signal for each 
location, the multiple receiver locations being spaced 
apart from one another and from the source transducer 
location. For each receiver location, the complex values 
of the received sonic signal are detected relative to a 
predetermined phase reference. The detected complex 
values of the received sonic signal may be processed to 
determine multipole sonic wave velocity and shear 
attenuation at discrete frequencies and to calculate ve- 
locity and attenuation as a function of frequency. 
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DISCRETE-FREQUENCY MULUPOLE SONIC 
LOGGING METHODS AND APPARATUS 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention is directed to logging of forma- 
tion shear in a borehole traversing an earth formation. 
More particularly, the present invention is directed to 
logging of formation shear using discrete frequency 
measurements, as v^ell as to logging in the frequency 
domain. 

2. Background Information 

It has long been an objective of the oil industry to 
obtain formation shear d ata fr om a sonic logging tool. 
See, for example, J. WHITE, THE HULA LOG: A 
PROPOSED ACX)USTIC TOOL, presented at the 
Society of Professional Well Log Analysts Conference, 
1967. Shear wave data has been found useful in lithol- 



10 



15 



the onset of the compressional and shear arrivals and 
the tube (Stoneley) arrival can be seen in these exam- 
ples, although the wave packets are not totally sepa- 
rated in time. The first arrival (the compressional wave) 
has traveled from the transmitter to the formation as a 
fluid pressure wave, has been refracted at the borehole 
wall, has traveled within the formation at the compres- 
sional wave velocity of the formation, and has traveled 
back to the receiver as a fluid pressure wave. The shear 
wave has traveled from the transmitter to the formation 
as a fluid pressure wave, has traveled within the forma- 
tion at the shear wave velocity of the formation, and has 
traveled back to the receiver as a fluid pressure wave. 
The mud wave (not strongly evident in the wavetrain 
examples of FIG. 1) has traveled directly from transmit- 
ter to receiver in the mud column at the compressional 
wave velocity of the borehole fluid. The Stoneley wave 
is of large amplitude and has traveled from transmitter 
to receiver with a velocity less than that of the compres- 



ogy and fluid indentification, porosity determination, M ^^j^ waves in the borehole fluid. 



measurement of rock elastic and inelastic properties, 
and as an adjunct to shear seismic data. See, for exam- 
ple, H. LESLIE and F, MONS. SONIC WAVEFORM 
ANALYSIS: APPLICATIONS, presented at the Soci- 
ety of Professional Well Log Analysts (SPWLA) An- 
nual Logging Symposium, New Orleans, La. June 
10-13, 1984. 

Traditional sonic borehole logging tools have been 
used to acquire formation shear data. Such tools employ 



25 



Various waveform processing techniques have been 
used to find and analyze some or all of the propagating 
waves in composite time-domain waveforms such as 
those shown in FIG. L For example, the slowness-time 
coherence technique (STC) uses a semblance algorithm 
to detect arrivals that are coherent across the array of 
receiver waveforms and to estimate their interval transit 
time. Applying this semblance algorithm to the wave- 



♦ '"^ "r"^- Tu i ,n forms of FIG. 1 produces the coherence map of FIG. 2, 

monopole transmitter which emits a broad-band 30. ft v j^.t. 



sound pulse, and one or more receivers which detect 
the pulse as it passes. 

A number of such sonic logging tools have been 
widely used, including Schlumberger's BHC* borehole 
compensated sonic tool, LSS* long-spaced sonic tool, 
and Array Sonic* digital sonic tool. (The symbol * 
indicates a mark of Schlumberger.) The waveform from 
each receiver of such a tool is recorded as a function of 
time to produce a sonic log in which the time required 



in which regions of large coherence correspond to the 
compressional, shear, and Stoneley arrivals. The apex 
of each region defines the slowness of that wave. This 
process is repeated for each set of waveforms acquired 
32 by the tool and is used to produce a log of slowness 
versus depth such as shown in FIG. 3. In a slow forma- 
tion (a formation having a shear velocity less than the 
borehole fluid velocity), the tool obtains real-time mea- 
surements of compressional, Stoneley and mud wave 



for a sound wave to traverse a given distance of the 40 velocities, but shear wave values must be derived from 
formation is plotted versus depth in the borehole. Such velocities, 
tools have not, however, proven entirely satisfactory Despite the progress made in obtaining formation 
for acquisition of formation shear data. For example, it shear data with the traditional monopole, time-domain 
is not always possible to obtain formation shear data in hogging tools, a number of difficulties remain. For ex- 
"slow" formations (formations having a shear velocity 45 ample, it has often been found that the amplitude of the 
less than the borehole fluid velocity). A review of how shear wave is insufficient for effective processing and 
shear data is obtained with these tools helps to under- analysis. And because the first shear wave arrival fol- 
stand why this is so. lows the compressional wave arrival, identifying the 

Propagation of the broad-band, monopole sotmd first shear arrival can be difficult or impossible, 
pulse in a borehole is governed by the mechanical prop-/5o" Accordingly, various sonic tools have been proposed 
ertics of several separate acoustic domains. These in( for dircctiy logging formation shear, particularly in 



elude the formation, the borehole fluid column, and the 
logging tool itself. Wide-band sound energy emanating 
from the monopole transmitter impinges on the bore- 
hole wall, establishing compressional and shear waves 
in the formation, surface waves along the borehole waU, 
and guided waves within the fluid column. The bore- 
hole wall, formation bedding, borehole rugosity, and 
fractures can all represent significant acoustic disconti- 



slow formations where the conventional tools having a 
monopole source cannot give shear wave logs directly. 
These tools employ dipole or other multipole (e.g., 
55 quadrupole) source and receiver transducer (rather 
than monopole transducers as in conventional sonic 
logging tools) for direct logging of shear in both fast 
and slow formations. 
A monopole source excites azimuthally symmetric 



nuities. The phenomena of wave refraction, reflection 60 acoustic waves around the borehole, as shown diagram- 



and conversion lead to the presence of many acoustic 
waves in the borehole when a sonic log is being run. 
Thus, many acoustic energy arrivals are seen by the 
receivers of a sonic logging tool. 



matically in FIGS. 4a and 4b. As a monopole wave 
propagates along the hole, the wave causes the borehole 
cross-section to bulge and contract symmetrically. In 
contrast, a dipole source can be synthesized from two 



The more usual acoustic energy arrivals at the recciv- 65 monopole sources of opposite polarity (Le., the pressure 



ers of a sonic logging tool having a wide-band mono- 
pole energy source are shown in the waveform exam- 
ples of FIG. 1. The distinct changes corresponding to 



of one source is positive when the other is negative and 
vice versa) placed on a plane perpendicular to the bore- 
hole axis. Since one side of a dipole source pushes the 
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fluid with a positive pressure and the other side pulls rection when the other segment contracts in the radial 

with a negative pressure, a dipole source acts as a point direction, and vice versa, thereby generating a dipole 

force moving the fluid transversely to the borehole wave in the borehole. The receiver transducers are of 

wall. Therefore, propagation of a dipole wave causes the same construction, acting as a pressure gradient 

the formation surrounding the borehole to flex side- 5 sensor at low frequencies. A voltage pulse is applied to 

ways, as shown diagrammatically in FIGS. 5a and 5b, the source transducer to produce a flexural wave in the 

A dipole source must move the borehole fluid to borehole, and the response of the receiver transducer is 

vibrate perpendicularly to the borehole axis. This can displayed in time domain. 

be achieved by moving a plate in the hole. One way to A fiirther proposal for low-frequency dipole shear 
move the plate is to use the concept of moving coUs in 10 logging is found in C KTTSUNEZAKI, A New Method 
a magnetic field, similar to loudspeakers. Another way for Shear-wave Loggings 45 Geophysics 10, at 1489-1506 
is to use piezoelectric bimorphic or monomorphic bend- (October, 1980). In that arrangement, a force applied to 
crs. A bimorphic bender is made of two thin piczoelec- a rigid body suspended in the borehole fluid produces 
trie plates of opposite polarity bonded together. As a pressure changes at front and rear surfaces of the body, 
voltage is applied across the plates, one plate extends 15 The borehole wall is excited indirectly by the fluid 
while the other contracts. Thus, the composite plate pressure changes. The clectrodynamic source trans- 
bends in response to applied voltage much as a bimetal- duccr employs a moving coil and a permanent magnet, 
lie plate does in response to temperature change. When As shown in U.S. Pat. No. 4,207,96 1 to Kitsunezaki, the 
supported at the edges, the bender will move the fluid moving coil drives a diaphragm which displaces the 
sideways as a dipole source. Other dipole source config- 20 fluid to excite dipole waves in the borehole. The recei v- 
nrations are also known, such as dual volume sources, ers arc geophone-like devices with a neutrally-buoyant 
expanding rods, rare earth transducers, split cylinders, body which moves relative to the logging tool in re- 
etc, sponse to passage of the asymmetnc dipole pressure 

A dipole receiver must be able to sense either pres- wave. A pulse is applied to the source transducer, and 

sure gradient, particle acceleration, particle velocity, or 25 the waveforms detected at the receivers are recorded in 

particle displacement. Measuring the differential output the time domain. The tool must be stationary while 

of two hydrophones wDl give the pressure gradient takmg measurements. An alternate transducer construe- 

Accelerometers and geophones will sense particle ac- tion proposed in U.S. Pat. No. 4,207,961 to Kitsunezaki 

celeration and velocity, respectively. Variable capaci- has diaphragms which move up and down such that a 

tance microphones can sense the displacement of the 30 pair of fluid chambers displace fluid in opposite direc- 

fluid particle vibrations. These transducers vary in their tions to generate a dipole wave, 

sensitivities and frequency characteristics. A modification of Kitsunezaki's clectrodynamic 

In contrast to the dipole transducers, which excite transducer, disclosed in U.S. Pat. No. 4,383,591 to 

and detect pressure gradient or particle vibrations, a Ogura, generates both monopole and dipole signals, 

quadrupole transducer can excite and detect waves 35 Since the original design by Kitsunezaki was imsuitable 

with 90-degree asymmetry. That is, as a quadrupole for logging shear in fast formations, a new design was 

wave travels up the borehole, the borehole vibration developed which uses a stationary coil and a moving 

will be like squeezing a paper cup in one direction. In magnet or a magnetic material block. The moving mag- 

the cross-sectional plane, the wall will squeeze in net is used as a hammer to drive either a rubber or a 
toward the axis in one direction and expand out away 40 metal plate. The frequency was extended to several kHz 

from the axis in the perpendicular direction. Thus, and the power output was increased by the new design, 

quadrupole transducers will excite and detect shear/- See K. TANAKA, S. INOUE and K. OGURA, DE- 

screw waves, in contrast to the shear/flexural of dipole VELOPMENT OF A SUSPENSION PS LOGGING 

transducers. SYSTEM»S SEISMIC SOURCE FOR HARD 

An early multipole logging tool concept was pro- 45 GROUND, presented at the 56th SEG Annual Meet- 
posed by J. E. White in 1967. See J. WHITE, THE ing, Houston, Texas, Nov. 2-6, 1986. 
HULA LOG: A PROPOSED ACOUSTIC TOOL, In 1984, results of time-domain dipole shear logging 
paper presented at the Society of Professional Well Log using piezoelectric bimorphic bender transducers were 
Analysts Conference, 1967, and U.S. Pat. No. 3,475,722 published. See J. ZEMANEK, F. ANGONA, D. WIL- 
to J. E. White, in th at pr oposal, a cluster^ 50 LIAMS and R. CALDWELL, CONTINUOUS 
trie transducers 1f^S!ea"^on eacn ofseveral"i?a3s ACOUSTIC SHEAR WAVE LOGGING, presented 
wmcn arc prised against the borehole wall. Tbejtrans^ at the 25th SPWLA Logging Symposium, New Or- 
ducers of each pad are oriented along mutuaU yor^og* leans, 1984. U.S. Pat. Nos. 4,516,228 to J. Zemanek and 



cnal axes and selectiv ely activated to aM^y35p^^|^ 4,649,525 to F, Angona ct al. disclose dipole transduc- 
borenole w'au m a' d^^df^i§n!nf ne looi "can 'fie 55 ers, including both circular and rectangular benders, 
operated in one of four different modes to produce For the circular bender, the circumference of the 
torsional, flexural, radial and axial vibrations, by suit- bender element is mounted to the sonde by means of a 
able application of voltage pulses to the source trans- rubber mounting ring. Two edges of the rectangular 
ducers. The flexural (sideways) mode of vibration is bender are moimted with rubber mounting strips. The 
that of a dipole source. 60 source and receiver transducers are of similar constnic- 

A subsequent dipole tool concept, which does not tion. The benders have a resonant frequency such that, 
require contact with the borehole wall, is disclosed in when a broad-band voltage pulse is applied to them, the r\ 
U.S. Pat. No. 3,593,255 to J, E. White. The transducers frequency of the resulting acoustic signal ranges from *V^^t^U.t.1^6^ 
in that arrangement are each formed of two half-cylin- about 1 kHz to 6 kHz with a predominant frequency of 
der segments of piezoelectric material (barium titanate), 65 about 3 kHz. Benders of lower resonant frequency are 
assembled to form a cylindrical body. The segments are also discussed. 

connected electrically such that the application of a Still other multipole logging tools have been pro- 
voltage causes one segment to expand in the radial di- posed. For example, U.S. Pat. No. 4,606,014 to Winbow 
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et al. discloses an acoustic dipole logging device in excited more than the flexural wave. Therefore, the 
which the transducers comprise cither a bimorphic dipole signals can be viewed approximately as a mixture 
bender clamped at one edge or a pair of bimorphic of shear and flexural wave, with a velocity dispersion 
benders clamped at opposing edges. The source trans- following that of the flexural wave. Since the drop in 
ducer is driven by a voltage pulse. The outputs of the 5 low-frequency amplitude of the shear wave is quite 
receiver transducers are band-pass filtered and pro- steep, it is difficult to obtain reasonable shear signal 
cessed in time domain. U.S. Pat. No. 4,649,526 to Win- strength without also exciting the flexural mode, 
bow et al. discloses an acoustic logging device with Second, to get formation shear, a dipole tool must be 
dipole or higher-order multipole transducers having operated at very low frequencies where the shear exci- 
multiple piezoelectric members of split-cylinder or bi- 10 tation can be very weak. 

morphic plate form. The source transducers are driven Third, road noise (the noise generated by the sonde 
by voltage pulses. traveling in the borehole) can be very significant at low 

U.S. Pat. No. 4,682,308 to Chung and U.K. Patent frequencies. Transverse vibration of the tool is gener- 
Application GB 2 158 581 A of Chen et al, filed May 2, ally detected by the dipole receivers. It has been found 
1985, disclose further multipole acoustic logging trans- 13 that road noise energy is mainly concentrated below 1 
ducers, including such arrangements as quadrupole kHz and that, in general, road noise increases with in- 
transducers having vertically-mounted magnetostric- creased logging speed. 

tive rods or piezoelectric stacks, or horizontally Finally, the hole conditions in soft rocks can be poor, 
mounted magnetostrictive rods or piezoelectric stacks frequency making the dipole signals very weak in soft 
which are attached to vibrating masses. 20 rock wells. Of all these difficulties, signal-to-noise ratio 

Several other transducer arrangements have also and wave dispersion arc of primary concern, 
been proposed, such as an eccentric weight vibrator and As noted above, sonic logging measurements have 
clamped geophone to measure shear in the 50 to 300 Hz traditionally been made in time domain: a broad-band 
frequency band (W. BENZING, EXPERIMENTA- sonic energy source excites propagation of sonic waves 
TION IN DOWNHOLE SHEAR WAVE VELOC- 25 in the borehole, and waveforms detected at an array of 
ITY MEASUREMENTS, presented at the SEG Annu- nsceiver locations spaced from the source in the bore- 
lal Meeting, Las Vegas, Nevada, Sept. 11-15, 1983), a ^ole are recorded as functions of time. Shear wave 
bender transducer having two stacks of rectangular velocity at the dominant frequency is then estimated by 
piezoelectric bars which are pre-stressed in a mounting waveform stacking techniques, such as semblance or 
assembly (G. NUNN and W. CURRIE, BENDER- 30 radon transforms. However, time-domain recording has 
BAR TRANSDUCERS FOR EARTH ACOUSTIC several disadvantages if the frequency domain results 
MEASUREMENTS presentation at the Interwell Scis- are desired. The waveforms must be sampled in small 
mic Surveying Workshop, Los Alamos, New Mexico, time steps over a long period, resulting in a large set of 
Mar. 24-26, 1988), and a dipole transducer using four data to transmit, store, and process. The truncation in 
rare-earth rods to excite four mass vibrators in the fre- 35 time can cause interference in the processing. In addi* 
qucncy range of 2 kHz to 4 kHz (S. COHIC and J. tion, the signal to noise ratio can be quite small for the 
BUTLER, Rare-^arth Iron Square Ring Dipole Trans- wide-band transient measurements. The noise problem 
ducer, 11 JOURNAL OF THE ACOUSTICAL SOCI- can be significant for low frequency data Oess than 1000 
ETY OF AMERICA Aug. 2, 1982). U.S. Pat. No. Hz) because ofthe noises generated by the sonde travel- 
4,700,803 to Mallett and Minear discloses a magneto- 40 ing in the borehole. 

strictive transducer similar to that in the aforemen- In accordance with the present invention, measure- 
tioned U.K. Patent Application GB 2 158 581 A of ments are taken in frequency domain, thereby avoiding 
Chen et al. the above-noted problems. A frequency domain mea- 

However, several difficulties are encountered in ob- suremcnt is desirable for a tool that gives a single domi- 
taining true formation shear from dipole waveforms. 45 nant arrival. This is particularly true if the arrival is 
First, it has been found that a dipole source excites not dispersive, as is the case with the shear/fiexural arrival 
only the de^ed shear waves, but also a borehole flex- of dipole waves. 

ural wave. A. KURKJIAN and S. CHANG, Acoustic In contrast to the traditional tiroc-domain logging 
Multipole Sources in Fluid-Filled Boreholes, 51 GEO- techniques, U.S. Pat. No. 3,330,375, issued July 11, 1967 
PHYSICS 1 (January. 1986), at 148-163, and U.S. Pat 50 to J.E. White proposes a form of acoustic well logging 
Nos. 4,698,792 and 4,703,460 to Kurkjian et al. The in which the propagation velocities of compressional, 
flexural wave is dispei^ve (its phase velocity varies shear, mud and casing waves are determined from the 
with frequency) and is excited at a greater amplitude expression velocity = frequency X wavelength by em- 
than the formation shear by a wide*band dipole source. ploying transmitter and/or receiver tuning techniques 
FIG. 6 illustrates the dispersive nature of the flexuraMS to determine the wavelength for a known frequency. A 
waves as a function of frequency. The flexural wave variety of such techniques arc disclosed, involving 
phase velocity approaches the formation shear velocity wavelength tuning by varying frequency and/or phase 
at low frequency and gradually becomes slower at shift. In all the techniques, frequency and/or phase shift 
higher frequencies. FIG. 7 shows one example of the is adjusted until an amplitude peak is observed, the 
relative excitation amplitudes of the shear wave and 60 frequencies which produce such amplitude peaks are 
flexural mode as functions of frequency. The direct noted, and wave propagation velocity is calculated 
formation shear is much weaker than the flexural wave from the expression given above. In one such technique, 
at higher frequencies, but stronger than the flexural the wavelength of the transmitted signal is tuned to a 
wave at low frequencies where the flexural wave phase fixed wavelength of the transmitting and receiving ar- 
velocity is cloisc to that of the shear velocity. (As shown 65 rays by varying its frequency. In another, the wave- 
in the example of FIG. 7, the flexural mode is more than length of the arrays is tuned to a fixed value of the 
20 dB stronger than the shear wave above a certain transmitted wavelength by varying the phase shift be- 
frequency.) At low frequency, the formation shear is tween adjacent transducer elements. In another, the 
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frequency of the transmitted signal or the phase shift 
between adjacent receiving transducers, or both, arc 
varied. 

The logging techniques of U.S. Pat. No. 3,330,375 
>vou]d have a number of disadvantages for shear/flexu- 
ral or shear/screw wave logging. For example, use of 
monopole sources significantly reduces the signal-to- 
noise ratio of the shear/flexuraJ waves, since much of 
the emitted sonic energy from the monopole source will 



8 



of the invention a multipole source transducer emits 
sonic energy at a location in a borehole of at least one 
discrete frequency which induces propagation of multi- 
pole sonic waves in the borehole; sonic energy of the 
multipole sonic waves is received at multiple multipole 
receiver locations in the borehole, the multiple receiver 
locations being spaced apart from one another and from 
the source transducer location; and for each receiver 
location, amplitude of the received sonic energy of the 
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propagatejn TOmprcssjonal^ rather 10 multipole sonic waves is detected and phase of the re 

^_ . ceived sonic energy of the multipole sonic waves rela- 

tive to a phase reference is detected. The source trans- 
ducer may be driven continuously at one or a plurality 
of discrete frequencies, or may be driven intermittently, 
or may be driven alternately from a plurality of sources 
which provide one or a plurality of discrete frequencies, 
in one preferred form, the source transducer is driven 
from a repetitive signal pattern to emit sonic energy at 
a plurality of discrete frequencies. The phase reference 
is preferably derived from the same signal source which 
drives the source transducer. 

Shear logs are preferably obtained by sorting the 
dispersion results after processing the discrete fre- 
quency data. The detected amplitude and phase of the 
received sonic energy of the multipole sonic waves may 
be processed to determine wave phase velocity and 
wave attenuation at discrete frequencies and to calcu- 
late velocity and attenuation as a function of frequency. 
The multipole sonic wave data are useful, for example, 
in lithology and fluid identification, porosity determina- 
tion, measurement of rock elastic and inelastic proper- 
ties, and as an adjunct to shear seismic data. A number 
of different techniques may be used to process the multi- 
pole sonic wave data, such as three receiver analytical 
solutions, least square fittings, two receiver phase differ- 
ence and amplitude ratio, Prony*s method, convenrional 
beamforming, and the maximum likelihood method. 



mode. Further, the disclosed method and apparatus are 
directed to determination of wave velocity, and do not 
provide for logging of amplitude and phase so that both 
velocity and attenuation can be determined as functions 
of frequency. Because of the dispersive nature of the 
shear/flexural mode (wave velocity varies with fre- 
quency), both amplitude and phase (or real and imagi- 
nary parts) of shear wave energy are needed to under- 
stand the dispersion. Furthermore, the disclosed 
method requires downhole tmiing of the frequency 20 
tnd/or the wavelength while acquiring the data. The 
method is not directly applicable to obtain the fre- 
quency dependent characteristics of dispersive waves, 
such as the shear/flexural or shear/screw waves. 

In another approach, U.S. Pat. No. 4,419,748, issued 25 
Dec. 6, 1983 to R. W. Siegfried, II., proposes a continu- 
ous wave sonic logging method in which a continuous 
sine wave at a single frequency is emitted and received, 
and a spatial Fourier transform is performed over the 
receiver array. The resulting spatial frequency coropo- 30 
nents are then used to indicate the velocities of various 
sonic paths. The logging method of U.S. Pat. No. 
4,419,748 would have several disadvantages for multi- 
pole shear logging. For example, the method requires a 
large number of receivers in order to facilitate the spa- 35 
tial Fourier transforms. Further, the disclosed method 
records the instantaneous values of the received signal; 
that is, it is recorded in time domain. Therefore, the 
measurement is subject to noise interferences. No im- 
provement in the signal to noise ratio is realized by the 40 
proposed method. Furthermore, the measurement can 
only be done one frequency per logging run (due to the 
time domain recording). This would require numerous 
logging runs for the dispersive waves, for which the 
wave characteristics are functions of frequency. Since 45 
the logging time is a costly factor in wire line logging 
services, the method is not practical for logging disper- 
sive waves. 

It is an object of the present invention to provide 
methods and apparatus for borehole shear/flexural log- 50 
ging in which the aforementioned disadvantages of the 
prior art are avoided. 

It is a further object of the present invention to pro- 
vide methods and apparatus for borehole shear wave 
logging in which amplitude and phase (or real and 55 
imaginary parts) of shear/flexural energy propagated in 
the borehole and surrounding formation are detected 
for use in determining parameters such as shear wave 
phase velocity and shear wave attenuation as functions 
of frequency. 

These and other objects of the present invention will 
become apparent from the description which follows 
with reference to the accompanying drawing. 

SUMMARY OF THE INVENTION 

The present invention is directed to methods and 
apparatus for sonic shear logging of boreholes using 
discrete frequency measurements. In a preferred form 
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BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 shows waveform displays of typical acoustic 
energy arrivals at the receivers of a sonic logging tool 
having a wide-band monopole energy source; 

FIG. 2 shows a slowness-time coh^ence map result- 
ing from the application of a typical semblance algo- 
rithm to the waveforms of FIG. 1; 

FIG. 3 shows a log of slowness versus depth for the 
data of FIGS. 1 and 2; 

FIGS. 4a and 46 illustrate diagrammatically in cross- 
sectional elevation and plan views, respectively, the 
azimuthally-symmetric acoustic waves excited by a 
monopole source in the formation surrounding a bore- 
hole; 

FIGS. 5a and $b illustrate diagrammatically in cross- 
sectional elevation and plan views, respectively, the 
azimuthally*asymmetric acoustic waves excited by a 
dipole source in the formation surrounding a borehole; . 

FIG. 6 is an example of a plot of the ratio of flexural 
mode phase velocity to shear wave velocity as a func- 
tion of frequency, illustrating the dispersive nature of 
flexural waves, and illustrating that the f lexural wave 
phase velocity approaches the formation shear velocity 
at^ low frequenc y and grad ually becomes slower-^ 
^her trequeiig lK! 

MG. 7 Is an example of a plot of flexural mode and 
shear wave amplitude as a function of frequency; 

FIG. 8 shows in schematic block form a first pre- 
ferred embodiment of a discrete-frequency dipole log- 
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ging apparatus for shear/flexura! logging in accordance in amplifiers 32 is passed to conventional telemetry 
with the invention; circuitry 50 for transmission to up-hole system 52 where 

FIG. 9 shows in schematic block form a modification the data is recorded and/or processed. Oscillators 28, 
of the embodiment of FIG. 8 which permits multiple, 30, 32, 34, mixer and 36, power amplifier 38, band-pass 
discrete frequencies to be emitted sequentially in bursts 5 filters 40, lock-in amplifier and phase-detector circuits 
for shear/flexural logging in accordance with the in- 42, 44, 46, 48 and telemetry circuitry 50 are preferably 
vention; built into sonde 10. Telemetry circuitry 50 preferably 

FIG. 10 shows amplitude, of a periodic signal as a transmits data output from lock*in amplifiers 36 to up- 
fimction of time, the signal having a discrete frequency hole recording and processing system 52 via cable 12 or 
spectrum as shown in FIG. 11; 10 any other suitable means. Alternatively, the data may be 

FIG. 11 shows amplitude of the discrete frequencies stored in the sonde and retrieved once the sonde has 
of the signal of FIG. 10; and been removed from the borehole. 

FIG. 12 shows in schematic block form a further Sourcetransducer 16maybeanytypeofdipole wave 
modification of the embodiment of FIG. 8 which per- source suitable for shear/flexural wave logging. Trans- 
mits multiple, discrete frequencies to be emitted simulta- 15 ducer 16 may be made ofpiezoelectric ceramics made in 
neously for shear/flexural logging in accordance with benders or cylindrical sections, magnctostrictive trans- 
the invention, ducers, or elcctrodynamic vibrators, or may be of any 

DETAILED DESCRIFTION OF THE other suitable construction. The preferred source trans- 

PREFERRED EMBODIMENTS moving coU elcctrodynamic transducers. 

2^ piezoelectric bimorph benders, or piezoelectric mono- 
FIG. 8 shows in schematic block form a first pre- morph benders. Source transducer 16 may be driven in 
ferred embodiment of a logging apparatus in accor- any manner suitable for exciting propagation of shear A 
dance with the present invention, useful for obtaining flcxural waves in the formation surrounding the bore- 
shear/flexural data at discrete frequencies. Referring to hole at one or more discrete frequencies, only some of 
FIG. 8, acoustic logging sonde 10 is shown suspended 25 ^iijch will be described herein, 
cable 12 in borehole 14 in the earth. Borehole 14 i^. 

Receiving transducers 20-26 are likewise any type of 
filled wi^ drmms mud or other fluids (not iHpstrated). jipolc transducers which detect pressure gradients or 
Loggmg sonde 10 preferably mcludcs dipole transmit- particle vibrations (acceleration, velocity or displace- 
ting, or source transducer 16 and a plurality of dipole ^^^y ^he preferred receiver transducers include dual 
receiving transducer 20-26 arrayed at locations spaced 30 hydrophones, benders, and elcctrodynamic transducers, 
apart from one another andfrom source transducer 16 ^hc source and reciving transducers should satisfy the 
along sonde 10. While FIG. 8 shows four reccivmg ^ specifications on temperature, pres- 

transducers, more or fewer uansducers may be used, maximum diameter, shock and vibration, and me- 

although It will be seen from the data proccssmg tech- chanical strength 

niques describe! below that at least two receiving uans- 35 frequency at which the flexural waves have near 

ducers arc preferred. velocity can be estimated by the foUowing for- 

Narrow-band transducer dnvers (such as discrete muJa- 
frequency oscillators) 28, 30, 32 and 34, generate sinu- 
soidal signals at one or more discrete frequencies fi, f2, 
f3> • * • ffl- If multiple discrete frequencies are produced, 40 fs: — ^ 
the signals are then mixed in signal mixer 36. The output 

signal from mixer 36 is then passed through amplifier 38 , . « ^ . , , . . , 

to drive source transducer 16. The signals may be con- ^^^^^^^ ^^^^ formation shear velocity and a is the 
tinuous. or may be transmitted in brief bursts of sufli- borehole radius. It is preferred tiiat the excitation fre- 
cient duration (for example, on the order of 10 or 20 45 f^'^^X ^ ^«Pt ?^ ^^w as possible to ensure that the 
cycles per burst) to approximate a continuous wave for fommtion shear is measurable. As a guidelme, the cir- 
purposes of exciting propagation of shcar/flcxural cumference of the borehole should be larger than a 
waves in the borehole at the desired frequency or fre- ^hear wavelength and smaller than twice the shear 
qucncies. When so driven, source transducer 16 induces wavelength. That is. the preferred frequency band is: 
propagation of shear/flexural waves in the earth forma- SO 
tion surrounding the borehole at the prescribed fre- Vs 2K, 

quency or frequencies. The shear/flexural waves are 2va <^< zna 

detected at receiving transducers 20-26. The signal 

output from each of receiving transducers 20-26 is scpa- where f is the frequency. V, is the formation shear ve- 
rated into a respective signal at each of the discrete 55 locity, and a is the borehole radius. For example, the 
frequencies of interest (fj, f2, fa. . . . fn) by a plurality of frequency range is about 1 kHz to 2 kHz for an 8-inch 
band-pass Alters (shown generally at 40). The output diameter hole with a 450 ^/ft shear slowness. The 
signal from each of band-pass filters 40 is supplied to a frequency range would be reduced to between 500 Hz 
respective one of lock-in amplifier and phase-detector and 1 kHz for a 16-inch diameter hole with the same 
circuits 42, 44, 46, 48, which detects the amplitude and 60 shear slowness. 

phase (or real and imaginary parts) of the received To summarize, the frequency for dipole shear logging 
waves at a respective one of the selected discrete fre- should be lower for larger holes and for formations with 
quencies fi, fs, fa. . . . f^. Each of lock-in amplifier and slower shear. A typical frequency band is from 500 Hz 
phase-detector circuits 42, 44, 46, 48. also receives a to 5 kHz, for shear slowness up to 450 /is/ft and bore- 
phase reference signal from a respective one of oscilla- 65 hole diameter up to 16 inches. Another consideration is 
tors 28, 30, 32, 34. The phase reference signals are pref- the shale alteration, where the formation near the weU is 
erably the signals generated by oscillators 28, 30, 32, 34 substantially softer than the unaltered formation. In this 
for driving source tranducer 16. Data detected by lock- case, the effective size of the borehole is bigger and the 
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frequency range employed should be lower than that purposes of the present invention, the Fourier spectrum 

used for a simple borehole in a faster formation. of the signal preferably has much higher energy in its 

The foregoing open-hole frequency considerations lower frequency components than in its higher fre- 

can also be applied to cased holes. In cased holes, it is quency components. This is because the upper fre- 

preferred to limit the frequency to a frequency lower 5 quency boundary for the energy of shear/flexural 

than that which gives a casing shear wavelength equal waves propagating in a borehole is about 5 KHz for 

to the casing circumference. And, if the casing is un- boreholesof a diameter typically encountered in oilfield 

bonded by the cement, the preferred frequency f should well logging. (It will be recalled that this upper frc- 

also be selected such that f< [(velocity of shear in the quency boundary is dependent on borehole diameter 

steel casing)/(2irxinean casing radius)). Tliis is a fre- 10 and the formation surrounding the borehole.) Greater 

quency at which the mean circumference of the casing energy at the low-frequency end of the spectrum will 

k one wavelength of the shear in the casing steel. thus provide greater signal-to-noisc ratio of the shear/- 

Bccause of road noise (the noise generated by the flexural waves. Selecting minimum and maximum fre- 

sonde traveling in the borehole), a high-power dipole qucncies of f„/„ and f««, respectively, then the perio- 

lourcc transducer is desirable for low-frequency opera- 15 dicity of the signal may be chosen such that 
tions (particularly below 1 kHz) in slow formation and- 

/or large hole environments. Stronger transducer out- t= if every u tn integer multiple of the 

put gives larger signal-to-noise ratio, and allows for mimmum frequency. 
faster logging speed. 

Other considerations in designing and operating a 20 ^ low-pass filter can be applied to the periodic signal 

sonde in accordance with the invention include sonde to eliminate the harmonics above fmax- Other (e.g., 

housing wave coupling, impedance matching, cffi- square-wave, saw-tooth, specially-designed pulse shape, 

ciency, sensitivity of receivers, etc. One way to reduce etc.) periodic signals may also be used for exciting dis- 

the sonde coupling is to reduce the mass of the vibrating crete frequency signals. 

clement of the source transducer, so that the mass of the 25 FIG. 12 shows in schematic block form a modifica- 

vibrating element is small relative to that of the station- tion of the embodiment of FIG, 8 which permits multi- 

ary sonde structure. It is also preferred to mechanically plc, discrete fr^uencies to be emitted simultaneously 

isolate the transducers from the rigid sonde housing for shear-wave logging in accordance with the inven- 

structure, such as by means of an isolation section. tion. Referring to FIG. 12, the output signal from wave- 

FIG. 9 shows in schematic block form a modification 30 form generator 70 is supplied to low-pass filter 72 which 
of the embodiment of FIG. 8 which permits multiple, passes frequencies f=fma;r. The filtered signal is then 
discrete frequencies to be emitted in sequential bursts passed through power amplifier 74 to drive source 
for shear/flexural logging in accordance with the in- transducer 16. The emitted sonic energy from source 
vention. A plurality of oscillators, for example three transducer 16 excites shear/flexural wave propagation 
oscillators 60, 62 and 64, generate sinusoidal signals at 35 in the borehole and the shear/flexural waves are re- 
respective frequencies suitable for exciting shear/flexu- ceived at each of the spaced-apart receiving transduc- 
ral waves m the borehole, for example at frequencies of ers, only one of which is shown (at 20) in FIG. 12. The 
I KHz, 3 KHz and 5 KHz. Suitable switching means 68 output signal from each of the receiving transducers is 
alternately passes multi-cycle bursts of the respective passed through a respective band pass filter (one of 
frequencies to power amplifier 38 to drive transducer 16 40 which is shown at 76 in FIG. 12) having a pass-band 
and to lock-in amplifiers 42-48 (FIG. 8) for use as phase frequency characteristic fmm<f<^mxu- The filtered 
references. The multi-cycle bursts arc of sufficient signal is then supplied to the input of a phase-sensitive 
length (for example, 10 to 20 cycles) to effectively simu- detector circuit (one of which is shown at 78 in FIG. 12) 
late a continuous signal for purposes of exciting shear/- which detects the phase of each frequency component 
flexural propagation in the borehole at the selected 45 of the filtered signal relative to a respective phase refer- 
frequencies. cnce signal, as well as the amplitude of each frequency 

A further modification of the embodiment of FIG. 8 component. The phase reference signals may be ob- 
may be made which permits driving the source trans- tained from the circuitry which drives source trans- 
ducer or transducers at a desired set of frequencies by ducer 16. Phase-sensitive detector 78 may be of any 
using a repetitive signal sequence. If the period of the 50 suitable analog, digital or hybrid construction. The 
sequence is r, then a discrete frequency spectrum at phase and amplitude data for each frequency is passed 
f==n/r is obtained, where n^ 1, 2, 3, etc. The amplitude by suitable means (for example, via telemetry to an 
distribution for the frequencies of interest can be de- up-hole site) from phase-sensitive detector 78 to a pro- 
signed by the shape of the waveform of the sequence. In cesser 80 which makes use of the data for any of a 
other words, the shape of the waveform can be de- 55 number of purposes, some of which are described be- 
signed such that a desired discrete frequency spectrum low. 

is obtained. This may be achieved by first specifymg the AS noted above, the upper frequency boundary for - 

amplitude distribution for the discrete frequencies, and the energy of shear/flexural waves propagating in a ?' 

then by summing the Fourier series of all the discrete borehole is about 5 KHz for boreholes of a diameter 

frequency terms with the pre-assigncd amplitude distri- 60 typically encountered in oilfield well logging. This 

, , ^PP«f frequency boundary is dependent on borehole 

One example of such a repetitive signal sequence is diameter and the formation surrounding the borehole, 

given in FIG. 10, which shows the amplitude of a peri- Thus, the discrete frequencies at which shear/flexural 

odic signal as a function of time. This periodic signal has waves are to be excited in the borehole preferably in- 

a discrete frequency spectrum at f„=0.5, 1. 2, 3, 4, 5 65 eludes a number of discrete frequencies up to the upper 

kflohcrtz. The Fouri er amplitude j^pY ear^^fre/^jipfipy ^ frequency boundary, for exainple 0 5 KHz 1 KHz, 2 

roversely proportional^to^f^ri'^us. the Fnurier sp^rfnim KHz, 3 KHz, 4 KHz and 5 KHz. Spacing of the rccciv- 

df the signal of FIG. 10 is as shown in FIG. 11. For in g transducers is not critical hnt ifp| -p ff?afa!v sucli t^^SL 
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the length of the receivine transducer array is not sub - or the real and imaginary parts. Since the optimal fre- 

/ stantiaBy less man the wavelen gth in order to accu - quency for the multipole tool can vary with formation 

fiitelydeterminetfae sheai/flaural wave phase veloc- shear and borehole diameter, it is preferred that wave- 

;ty- . forms are recorded at more than one frequency. Prefer- 

Further modification of the embodiments of FIGS. 8, 5 ably, the low and high frequency limits are first se- 
9 and 12 is also possible within the spirit and scope of lected. For example, these may be 100 Hz and 5000 Hz 

the present invention. For example, quadrupole or for a dipole tool. Several discrete frequencies are se- 

other higherorder multipole transducers may be used. lected between the frequency limits. The data at each 

Several differences between dipole and quadrupole frequency is preferably then processed to obtain shear 

shear logging should be noted. Dipole logging uses lo velocity using one ofthe methods described below. By 
dipole transducers which cxdte and record pressure' sorting the frequency-dependent velocity and ampli- 
gradient orpirtide vibration. Quadnapole logging uses\Ltude data at the selected Muendes. the best choice of 

qtiadrupole transducers which excite and detect waves formation shear parameters can then be made by cross- 

with 90slepee asymmetry. A5 a quadrupole wave trav- checking with the theoretical models, 

els up the borehole, the borehole vibration will be hke is por example, one approach to obtain shear velocity is 

«!ueezing a paper cup in one direction. In the cross^- ^ select the velocity and ampUtude at the lowest possi- 

tiona^ plane, the wall wiU squeeze m toward the axis m frequency which has a reasonable signal-to-noise 

one direction and expand out away from the axis m the Qne way to judge the signal-to-noi^tio is to 

perpendicular direction. The quadrupole receivers are ^ ^ ^, frequencies, 

d^ed to sense this kind of vibration. - 20 since the velocity as a function of frequency is smooth 

_ Second, the quadrupole transducers will exate and „j predictable by theory, the quality of the signal, or 

detect shear/sCTew w.ves^ rather than the shear/flexu- ^le signal-to-noise ratio is sm^, when the s^ess 

ral waves of dipole transducers, Shear/flexural wave, JtL „r a.™ fi!^ ♦^^j „4 

J 1.- t. J 1 estimates are out of line from the trend produced at 

t^J^^ltZ n 1 »f8^^f^-«.'<*«^ °"»^?P° other frequencies. Another way is to check the coher- 

niird. th^ frequency selection for quadrupole log- °^ ''?!L*^""v ^ 

ging be dilfwent from that for dipole logK For ^ ^'7^%^ ' °" n n"!? 

iuadnipolc logging, the rule of thumbwill Kt the ffvSJ^''r°°T''i^^ f '''T*'^ 

drcuXenceTuie borehole should be Urger than 30 J^S^^tl l'^T °P^ " 

two shear wavdengths and smaller than thri shear ^P^n^L"' , rwAVr l^lr^' 

wavelengths. That fc, the preferred frequency band is: ^"^"^I! -fL^^ » ^ ^"Ao«^tT ' 

-1 J Sources m Flutd-Filkd Boreholes, 5\ GEOPHYSICS I, 

January, 1986 incorporated herein for all purposes by 
2H. this reference), at 148-163, or from FIG. 6, the low-fre- 

Zwo 2iro 35 quency limit of the shear/flexural or shear/screw veloc- 

... ^, , ■ .1. r ity is the formation shear velocity. Another approach is 

where f is the frequency V, is the formation shear ve- ^ . of fiequency 

lo«^y, and a B the borehole radius from the borehole measurements, and then search for 

Measurmg multipole some waves (e.g.. dipole shear/- j^^^^^j^^ ^^ear vdodty by using curve-fitting 

flexural waves or quadrupole shear/screw waves) m.40 methods, such as a least-square method or a generalized 

linear-inversion method. The theory has beel explidtly 
pie, measurement >« frequency domam give higher 5^ ^ KURKJIAN and S. CHANG, Acoustic 

l^t^^ ^^. ti^l^^^^J^^^^ '^''Mpole Sources in Fluid-Filled Boreholes, 51 GEO- 
ments. This is because the signal-to-noise ratio is m- nirvcu-'c i i.....,^. losc ias ti,_ij~.i., ,v:- 
versely proportional to the bandwidth ofthe signal, and 43 PHYSICS 1, January 1986^t 148-163. Exphcitly, th.s 
the bandwidth of the signal can be made ^reinely P^^^"f' is expres«d m two ways; 
narrow when workmg with discrete frequendes. Also. "'"^ '"'^'^^^ v^^such that: 

as there is no need to record waveforms, the quantity of 

data to be transferred up-bolc via the telemetry and j ^^^^^ _ vyf(/3|2i,nmuinum, 

cable is much less than that for tranditional time-domain 50 ^="1 ^ 

sonic waveforms. Thus, the logging speed can be in- 
creased. Additionally, as there is no need to transform where vpis the measured flexural- or screw-wave or 
the waveforms from time domain to frequency domain, other higher-order multipole wave velocity, vjf is 
the processing is simplified. Further, transducer calibra- the theoretically-calculated flexural- or screw-wave 
tion at specific frequencies is more reliable than is trans- 53 or other higher-order multipole wave velocity, and 
ducer calibration for the broad-band time domain is a weighting function which can be chosen to be 

sondes. Measurement in the frequency domain, as well one or a function of wave ampUtude at the i-th fre- 
ts the use of discrete frequencies, are believed to make quency. 

the methods and apparatus of the present invention 2. Find the shear velocity Vxj^carSuch that: 
more accurate than traditional time-domain sonic log- 60 
ging for obtaining the frequency-dependent information 

of multipole sonic (c,g., shcar/flcxural or shear/screw) ^ . ^^^^ minimum, 

waves. 1=1 N 

The data obtained with the methods and apparatus of 
the present invention are the complex values of the 65 The approach to solving these problems need not be 
dipole (or other multipole) signal at each receiver for described here, since it is well known in general. See, 
each frequency. It is well known that a complex value for example, R. HAMMING, NUMERICAL METH- 
can be expressed in two forms: the ampUtude and phase, ODS FOR SCIENTISTS AND ENGINEERS, 2d 
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Ed., 1973 (incorporated herein for all purposes by this plex value) of the n-th receiver located at 2„, and if the 

reference), at pages 68-70 (Newton's Method), and R. receivers are equally spaced at interval of Az, then 

WIGGINS, The General Linear Inverse Problem: Impli- Equation (3a) can be rewritten as 
cations of Surface Wa\e and Free Oscillations for Earth 

Structure, 10 REVIEWS OF GEOPHYSICS AND 5 Pn^A e-*"- (4) 
SPACE PHYSICS 1, 1972 (incorporated herein for all 

purposes by this reference), at pages 251-285 (general Taking the data from only three receivers, and defm- 

linwr mversc) . ^ ^ ing A'=A e^*^, B'-B t'^^ and G-e'*^, then th6 

Those of skill m the art will recognize that the mini- foUowing equations result: 

mization can be done over other norms of the differ- iq 

ences between the measured and theoretical velocity Px-A'G-^-^bg (5a) 
functions (the norms used above are square and linear 

norms) with proper choices of the weighting ftinctions. Z^^^'+J* (5b) 

The measurement can be made simultaneously for all 
the frequency points. The frequency oscillators can be u Py^A'G-k-BG-^ (3c) 
mixed and employed to drive the transmitter simulta- 
neously at all the desired frequencies. Narrow-band Eliminating A' and B', Equations (5a)-(5c) can be 
filters arc applied to the received signals to separate the reduced to a quadratic equation of G. The solution of 
data at the different frequencies. the quadratic equation is 

In the following, selected prior art processing tech- ^ 

niqucs are described which may be used to extract infor- _____ m 

mation from the measured data. Other processing tech- -k- H>± N(/>, + Pii^ -aPi^ 

niques may also be employed. ^ 2^ 

L TWO RECEIVER TECHNIQUE 
. . 1 . w v 1 . J. . 25 Once G is solved, k can be obtained from the cxpres- 

As a wave travels m the borehole m one dir«mon the ^y^^^^y ^he wave phase velocity, v. and 

detected-signal function can be expressed m the follow- ^3^, attenuation constant, a, can then be calculated 

™^ °"°* from Equation (3b). Solving for A' and B' from Equa- 

p^zS-A(/)e-^€^^(f^^^ (I) following expressions are 

30 obtained: 

where P^ is the detected signal (representing pressure 

gradient or other received quantities such as particle Pj - PjC (7a) 

displacement, velocity, acceleration, etc.); A is the am- " 1 - 

plitude function which is independent of transmitter- 
receiver spacings; z is the transmitter-receiver spacing; B^Pi^A' (7b) 
f is the frequency; v is the wave velocity; and a is the 

attenuation constant. The wave reflection coefficient can be obtained from 

If the complex detected-signal function is measured at the ratio of B' and A' with proper exponential factors, 
two receiver locations, zi and Z2, the wave velocity v(0 

and attenuation constant a may be calculated according ^ 3' ARRAY TECHNIQUES 

to the following formulas: Any of a number of known techniques may be used to 

urt-i-r h«*fPj n\ h*^pj process data from an array of receivers, such as: least 

A)]n'^^z2) (2a) method, Prony's method, maximum likelihood 

method, and others. 

ii^\a[<MPdnJ^^P^^iJ>M^i-^ ab) In the least square method, A and B and k are found 

such that 

The two receiver method has one constraint; the 
wave must be traveling in only one direction. If there ^^^y 

arc significant reflection signals from tool housing ^'^ 1^ {abs{Pa - AtP"» ~ Sc-^**)]^ 

structure, bed boundaries, or fractures, then other pro- "~ 

cessing procedures should be used which can include . . , „ * * 1 i_ * • 

both forward and backward traveling waves. minimizai JHere. the total number of receivers m the 

array is N. The least square method is feasible if the 
2. THREE RECEIVER TECHNIQUE . waveforms have only one arrival at the receivers (so 

Waves can be reflected, cither by the borehole/for- " Equation (4) is valid). If there are other arrivals, 
mation inhomogcneity or by the sonde housing struc- ^ compressional, fluid, tool arrivals and other 

ture. With the shear wave traveling up and down the waves, then special spectral estimation techniques (Pro- 
borehole, the detected-signal function can be expressed ^y'^, method, conventional beam-forming, maximum 
as; likelihood method, etc.) are preferred, The wave phase 

^ velocity v may be calculated frbtn' tfie"^^ 

v==(2vy)/Real Part (A), 

where k is the complex wave-number of the wave. 

and the wave anenuation constant a may be calculated 
Mnif/^Hfa (3b) 65 from the expression: 

If Pfl is deflned as the detected signal, representing imaginary Pux (k). 

pressure gradient or other received quantities (a com- 
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In applying Prony's method, the receiver response at 
a faed frequency b first assumed to be 



(9) 



A linear predictor is defined as h(z) so as to satisfy the 
relationship 



where the notation, means convolution in z. If H(k) 
is the Fourier transfonn of h(z) in the wave*numbcr 
domain. Equations (9) and (10) can be satisfied only 
when H(k) is zero at k^km- 

In Prony's method, the spatial variables V/z) and h(z) 
are first truncated to finite intervals of z. Ilic linear 
equations of Equation (10) are then used to solve for 
h(z) (assuming h(0) 1. The function h(z) is transformed 



sources are fired sequentially so that one estimate is 
from the upward-going waves and the other is from the 
downward-going waves. These two estimates (from the 
upward- and downward-going waves) are averaged to 
obtain the borehole-compensated results. Another way 
to achieve borehole compensation is to use depth- 
derived borehole compensation with only one source 
transducer. This is done by a synthesized array from a 
sequence of tool locations by picking the appropriate 
receiver so that the receiver location of the synthesized 
array is stationary. By estimating the formation quanti- 
ties from the synthesized array and from the original 
array data, and by averaging the two estimates at a 
proper depth location, a borehole-compensated result 
may be obtained. 

A further application is to derive the shear velocity 
profile as a function of radius from the borehole axis in 
the zone near the borehole frxim the collection of the 



13 



to H(k). Solving the k's for the equation H(k)=0, the 20 dispersion data at discrete frequencies. This can be done 

estimates of km*s are obtained. The wave phase velocity in the following way. First, the flexural wave velocity 

and attenuation are then obtained from the km's accord- versus frequency is estimated from the multipole log- 

ing to Equation (3b). The amplitude functions Am*s can ging tool signals and the processing described herein, 

also be estimated by least square fitting, (This is differ- Second, a theoretical model of acoustic waves in a bore- 

ent from the least square method described above. Here, 25 hole with radial velocity profile is used to calculate the 
the calculated values of km are used in the fitting.) Fur- 
ther information on Prony*s method may be found, for 
example, at S. L. MARPLE, JR., DIGITAL SPEC- 
TRAL ANALYSIS, Prentice Hall, 1987 (incorporated 

herein for all purposes by this reference), at pages ij^"*^^"*!*^ w j w ^ *u ^ v 

303-349 ' *- o 3U The search can be done by curve-fittmg methods such 

as a least square method or a generalized linear inver- 



flexural wave profile (for example, as described below). 
Then, one searches for the velocity profile until the 
difference between the measured and calculated (from 
theory) flexural velocity versus frequency is minimized. 



In the conventional beamforming method, the func- 
tion p(k) is calculated from the expression: 



(11) 



35 



where R=[P]p?]* is the covariance matrix from the 
measured data, and 



(H) 
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sion method. The theory can be derived from A. 
KURKHAN and S. CHANG, Acoustic Multipole 
Sources in Huid-Filled Boreholes, 51 GEOPHYSICS 
Jan. 1, 1986 incorporated herein for all purposes by this 
reference), at pages 148-163, or from Appendix A of L. 
BAKER and G. WINBOW, Multipole P-wave Logging 
in Formations Altered by Drilling, 53 GEOPHYSICS, 
September, 1988 (incorporated herein for all purposes 
by this reference), at pages 1207-1218. Explicitly, this 
procedure is expressed in two ways: 
1. Find the shear velocity as a function of radius v/r) 
such that: 



45 

The superscript ♦ means complex conjugate and 
transpose of the matrix or vector. The values of the 
wavenumbers km are obtained by searching the loca- 
tions of the maxima of p(k). 

In the maximum likelihood method, the function L(k) ^ 
is first calculated as 



flfcr[G«(A)iJ-lG(*)p 



(12) 



55 



where G(k) and R are in the form described above. The 
values of k^ are obtained by searching for the locations 
of the local maxima of L(k). Further information on the 
maximum likelihood method may be found in Capon, 
Greenfield and Kolker, Multidimensional Maximum 60 
Likelihood Processing of a Large Aperture Seismic 
Array. 55 PROCEEDINGS OF THE IEEE, 1967 
(incoTporated herein for all purposes by this reference), 
at pages 192-21L 

Borehole compensation may be. included in the tool 65 
or the processing algorithms. One way to achieve bore- 
hole compensation is to use two source transducers, one 
above and the other below the receivers. The two 



/V 



where v/jis the measured flexural wave velocity, yjj* 
is the theoretically calculated flexural wave velocity, 
and w/is a weighting function which can be chosen to 
be one or a function of wave amplitude at the i-th 
frequency. 

2. Find the shear velocity as a function of radius v/r) 
such that: 



^ VI* ^1 i* niimmnni. 

/«1 N 

The approach to solving these problems need not be 
described here, since it is well known in general See, 
for example. R. HAMMING, NUMERICAL METH- 
ODS FOR SCIENTISTS AND ENGINEERS, 2d Ed. 
(1973), at pages 68-70, and R. WIGGINS, The General 
Linear Inverse Problem: Implications of Surface Wave and 
Free Oscillations for Earth Structure^ 10 REVIEWS OF 
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GEOPHYSICS AND SPACE PHYSICS 1 (1972), at 
pages 251-285 (general linear inverse). 

It is noted that while many of the spectral*analysis 
processing methods described above have in the past 
been employed in efforts to transform time-domain 5 
signals into frequency domain (e.g., in efforts to esti- 
mate frequency domain velocities from time-domain 
data), it believed that they have not heretofore been 
applied to the estimation of velocities from frequency- 
domain measurements acquired in accordance with the 10 
present invention. 

Those of skill in the art will recognize that discrete- 
frequency, multipolc-wave logging tools other than the 
dipole and quadrupole tools described herein can also 
be configurod within the spirit and scope of the present 15 
invention. 

The preferred embodiments described above are not 
intended to be limiting, but are instead intended as 
merely illustrative of the present invention. Those of 
skill in the art will recognize that many modifications 20 
may be made in the disclosed embodiments without 
departing from the spirit and scope of the present inven- 
tion as defined by the following claims. 

What is claimed is: 

1. A method of borehole logging, comprising the 25 
steps of: 

(a) emitting sonic energy from a multipole source at a 
first location in a borehole, said sonic energy being 
of at least one discrete frequency which induces 
propagation of multipole sonic waves in the wall of 30 
said borehole; 

(b) receiving sonic energy of said multipole sonic 
waves at multiple receiver locations in said bore- 
hole to produce a received sonic signal for each 
location, said multiple receiver locations being 35 
spaced apart from one another and from said first 
location; 

(c) for each said receiver location, detecting the com- 
plex values of the received sonic signal; and 

(d) processing said detected complex values of the 40 
received sonic signals to determine phase velocity 
of said multipole sonic waves in said borehole as a 
function of frequency. 

2. The method of claim 1, wherein step (a) comprises 
continuously emitting said sonic energy at said discrete 45 
frequency. 

3. The method of claim 1, wherein step (a) comprises 
intermittently emitting said sonic energy at said discrete 
frequency. 

4. The method of claim 1, wherein step (a) comprises 50 
simultaneously emitting said sonic energy at a plurality 
of discrete frequencies. 

5. The method of claim 1, wherein step (a) comprises 
repetitively emitting said sonic energy at a plurality of 
discrete frequencies below 6 kHz. 55 

6. The method of claim 1, further comprising the step 
of repeating steps (aHc) for each discrete frequency. 

7. The method of claim 1. further comprising the step 
of deriving a predetermined phase reference from said 
emitted sonic energy, and wherein said step of detecting 60 
the complex values of the received sonic signal com- 
prises the step of detecting said complex values of the 
received sonic signal in relation to said predetermined 
phase reference. 

8. The method of claim 1, further comprising the step 65 
of processing said detected complex values of the re- 
ceived sonic signals to determine attenuation of said 
multipole sonic waves in said borehole. 



20 

9. The method of claim 8, further comprising the step 
of processing said detected complex values of the re- 
ceived sonic signals to determine attenuation of said 
multipole sonic waves in said borehole. 

10. The method of claim 6, further comprising the 
step of processing said detected complex values of the 
received sonic signals to determine attenuation of said 
multipole sonic waves in said borehole as a function of 
frequency. 

11. The method of claim 1, wherein step (d) com- 
prises the step of calculating multipole sonic wave 
phase velocity v and multipole sonic wave attenuation 
constant a from said detected complex values at two of 
said receiver locations and from the expressions: 

y))l/2ir/2i-i2) 

where the shear-wave detected-signal function is ex- 
pressed as: 

and where Pgis the detected-signal function; A is ampli- 
tude as a function of frequency and is independent of 
spacing between said emitting location and said receiver 
locations; z is the spacing between said emitting loca- 
tion and a receiver location; f is the discrete frequency 
at which the multipole sonic wave energy is propa- 
gated; and zi and Z2 are respective said receiver loca- 
tions. 

12. The method of claim 1, wherein step (d) com- 
prises the step of calculating multipole sonic wave 
phase velocity v and multipole sonic wave attenuation 
constant a from said detected complex values at three 
of said receiver locations spaced at intervals Az and 
from the expression: 

A=2ir(rA')-l-/a=(ln G)/(/Ar) , 

where 

m ■ — 

and where k is the complex wave-number of the re- 
ceived multipole sonic wave, f is the frequency of the 
received multipole sonic wave, and is the complex 
detected signal of the receiver located at an n-th one of 
said receiving locations, 

13. The method of claim 12, further comprising the 
step of calculating a multipole sonic wave reflection 
coefficient from the ratio of B' and A' in which 

Pz - P3C 

^ = - A'. 

14. The method of claim 1, further comprising the 
step of determining values for A and B and k such that 
the expression: 
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is miiumixed and where N is a plural number of said 
receiver locations, k is the complex wave-number of the 
received multipole sonic wave, Zn is an n-th one of said 
receiving locations, and is the complex detected 
signal of the receiver located at Zn and is expressed as: 

15. The method of claim 14, further comprising the 
step of calculating the multipole sonic wave reflection 
cocfQcient from the ratio of said determined values for 
A and B. 

16. The method of claim 14, wherein step (d) com- 
prises the step of calculating the multipole sonic wave 
phase velocity v from the expression: 

v=(2ir/VRaI Part (Jt). 



21. The method of claim 1, further comprising the 
steps of 

(a) calculating the function L(k) as 



10 
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where 



^krN 



17. The method of claim 14, further comprising the 
step of calculating the multipole sonic wave attenuation ^5 
constant a from the expression: 

a=Iiiugiiury Put (k). 

18. The method of claim 1, further comprising the 
steps of 

(a) defining the detected signal P(z) at a fixed fre- 
quency to be 



30 
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and where R=[P][P]* is the covariance matrix 
from the measured data, P» is the complex receiver 
response of the receiver located at Zn, N is the 
number of receivers in the array, and the super- 
script • denotes conjugate transpose of the matrix 
or tlie vector; and 

(b) determining values of ]Cm by searching for the 
locations of the local maxima of L(k). 

22. The method of claim 1, further comprising the 
steps of 

(a) calculating the function p(k)=G*(k) R G(k), 
where 



(b) defining a linear predictor h(z) so as to satisfy the 
relationship 



40 



and 



where the notation, (j), means convolution in z; 

(c) defining H(k) as the Fourier transform of h(z) in 
the wave-number domain, where said definitions of 45 
receiver response P(z) and said linear predictor 
h(z) can be satisfied only when H(k) is zero at 

(d) truncating spatial variables P(z) and h(z) to finite 
intervals of z; 50 

(e) solving the linear equations of the relationship 



W 1 A(z)==0 
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for h{z), assuming h(0)=l; 
(0 transforming the function h(z) to H(k); and 
(g) calculating values of k for the equation H(k) =0 to 

obtain estimated values for km- 
19. This method of claim 18, wherein step (d) com- 
prises the step of calculating the multipole sonic wave 60 
phase velocity and multipole sonic wave attenuation 
from said values for km from the relationship 



20, The method of claim 18. further comprising the 
step of estimating said amplitude functions Am by least 
square fitting. 



65 



and where R=[P][P]* is the covariance matrix 
from the measured data, Pr is the complex receiver 
response of the receiver located at Zn, N is the 
number of receivers in the array, and the super- 
script * denotes conjugate transpose of the matrix 
or the vectoi^ and 
(b) determining values of wavenimibers k^by search- 
ing for the locations of the maxima of L(k). 

23. The method of claim 1, wherein step (d) com- 
prises the step of deriving a radial shear velocity profile 
for a zone near said borehole. 

24. The method of claim 23, wherein said step of 
deriving a radial shear velocity profile comprises: 

(a) estimating a flexural wave velocity versus fre- 
quency profile in a zone surrounding said borehole 
from said detected complex values; 

(b) calculating a theoretical multipole wave velocity 
versus frequency and the shear velocity profile for 
said zone; and 
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(c) searching for a flexural wave velocity versus fre- 
quency profile for said zone by curve-fitting until 
the difference between said estimated profde and 
said calculated profile is minimized. 

25. The method of claim 1, wherein step (d) com- 
prises the step of determining shear velocity in a zone 
near said borehole by selecting the multipole sonic 
wave phase velocity at the lowest said discrete fre- 
quency for which a reasonable signal-to-noise ratio is 
obtained. 

26. The method of claim 1, wherein step (d) com- 
prises the step of deriving a shear velocity value for a 
zone near said borehole. 

27. The method of claim 26, wherein isaid step of 
deriving a shear velocity value comprises: 

(a) estimating a multipole sonic wave velocity versus 
frequency profile in a zone surrounding said bore- 
hole from said detected complex values; 

(b) calculating a theoretical multipole sonic wave 



10 
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velocity of said multipole sonic waves in said bore- 
hole as a function of frequency. 

31. The apparatus of claim 30, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
continuously emit sonic energy at said discrete fre- 
quency. 

32. The apparatus of daim 30, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
iBtermittently emit sonic energy at said discrete fre- 
quency* 

33. The apparatus of claim 30, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
simultaneously emit sonic energy at a plurality of dis- 
crete frequencies. 

34. The apparatus of claim 30, wherein said emitting 



velocity vereus frequency profUe for said zone; and 20 /^^P"?^ a multipole source transduder and 



(c) searching for a shear wave value for said zone by 
curve-fitting until the difference between said esti- 
mated profile and said calculated profile is mini- 
mized. 

28. The method of claim 27, wherein said step of 25 
searching for a shear wave value comprises the step of 
finding the shear velocity VjA^rSUch that: 



2 



N 
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where v/7 is the measured flexural- or screw-wave or 
other higher-order multipole wave velocity, vjf is the 
theoretically-calculated flexural- or screw-wave or 
other higher-order multipole wave velocity, and w/is a 35 
weighting function which can be chosen to be one or a 
function of wave amplitude at the i-th frequency. 

29, The method of claim 27, wherein said step of 
searching for a shear wave value comprises the step of 
finding the shear velocity v,Aeor such that: 40 



1=1.. 
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where v/r is the measured flexural- or screw-wave or 
other hi^er-order multipole wave velocity, vjt* is the 
theoretically-calculated flexural- or screw-wave or 
other higher-order multipole wave velocity, and w/is a 
weighting function which can be chosen to be one or a 50 
function of wave amplitude at the i-th frequency. 
30, Apparatus for borehole logging, comprising: 

(a) multipole source means for emitting sonic energy 
at a first location in a borehole, said sonic energy 
being of at least one discrete frequency which in- 55 
duces propagation of multipole sonic waves in the 
wall of said borehole; 

(b) multipole detector means for receiving sonic en- 
ergy of said multipole sonic waves at multiple re- 
ceiver locations in said borehole to produce a re- 60 
ceived sonic signal for each location, said multiple 
receiver locations being spaced apart from one 
another and from said first location; 

(c) means for detecting, for each said receiver loca- 
tion, the complex values of the received sonic sig- 65 
nal; and 

(d) means for processing said detected complex val- 
ues of the received sonic signals to determine phase 



means for driving said multipole source transducer with 
a repetitive signal such that said source transducer emits 
sonic energy at a plurality of discrete frequencies below 
6 KHz. 

35. The apparatus of claim 30, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
emit sonic energy in sequential bursts at a variety of 
discrete frequencies, and wherein, for each said bursty 
said multipole detector means detects for each said 
receiver location amplitude of the received sonic en- 
ergy of said multipole sonic waves and phase of the 
received sonic energy of said multipole sonic waves 
relative to a phase reference. 

36. The apparatus of claim 30, wherein said multipole 
detector means detects for each said receiver location 
amplitude of the received sonic signal of said multipole 
sonic waves and phase of the received sonic signal of 
said multipole sonic waves relative to a phase reference 
derived from said emitted sonic energy. 

37. The apparatus of claim 30. further comprising 
means for processing said detected complex values of 
the received sonic signals to determine attenuation of 
said multipole sonic waves in said borehole. 

38. The apparatus of claim 30, further comprising 
means for processing said detected complex values of 
the received sonic signals to determine attenuation of 
said multipole sonic waves in said borehole as a function 
of frequency. 

39. The apparatus of claim 30, wherein said process- 
ing means comprises means for deriving a radial shear 
velocity profile for a zone near said borehole. 

40. The apparatus of claim 30, wherein said process- 
ing means comprises means for deriving a value of shear 
velocity for a zone near said borehole. 

41. A method of borehole logging, comprising the 
steps of: 

(a) emitting sonic energy from a multipole source at a 
first location in a, borehole, said sonic energy being 
of at least one discrete frequency which induces 
propagation of multipole sonic waves in the wall of 
said borehole; 

(b) reccivbg sonic energy of said multipole sonic 
waves at multiple receiver locations in said bore- 
hole to produce a received sonic signal for each 
location, said multiple receiver locations being 
spaced apart from one another and from said first 
location; 
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(c) for each said receiver location, detecting the com- where 
plex values of the received sonic signal; and 

(d) processing said detected complex values of the ______ 

received sonic signals to determine attenuation of (/i + ± N (Pi + Pji^ - ap^ 

said multipole sonic waves in said borehole as a 3 

function of frequency. 

♦2. The method of clahn 41, wherein step (a) com- and where k is the complex wave-number of the re- 
prises continuously emitting said sonic energy at said ccived multipole sonic wave, f is the frequency of the 
discrete frequency. received multipole sonic wave, and P„ is the complex 

43. The method of dahn 41, wherein step (a) com- 10 detected signal of the receiver located at an n-th one of 
prises intermittently emitting said sonic energy at said said receiving locations. 

discrete frequency, 52. The method of claim 51. further comprising the 

44. The method of claim 41, wherein step (a) com- step of calculating a multipole sonic wave reflection 
prises amultaneously emitting said sonic energy at a coefficient from the ratio of B' and A' in which 
plurality of discrete frequencies. 15 

45. The method of claim 41, wherein step (a) com- ft _ p, g 
prises repetitively emitting said sonic energy at a plural- A' = . 
ity of discrete frequencies below 6 kHz. i - 

46. The method of claim 41, further comprising the ff A', 
step of repeating steps (aMc) for each discrete frc- 20 

quency. 53. jhc method of claim 41. wherein step (d) com- 

47. The method of claim 41, wherein step (d) further prises the step of determining values for A and B and k 
comprises the step of processing said detected complex j^ch that the expression: 

values of the received sonic signals to determine phase 
velocity of said multipole sonic waves in said borehole. 25 

48. The method of claim 41, further comprising the J? = "1 (a&j - At^ - se^f 
step of deriving a predetermined phase reference from ""^ 

said emitted sonic energy, and wherein said step of . . . ^ , ^ ^, . • , , ^ , , 

detecting the complex values of the received sonic sig- miBumred and where N is a plural number of said 
nal comprises the step of detecting said complex values 30 receiver locations, k is the complex wave-number of the 
of the received sonic signal in relation to said predeter- received multipole some wave. z„ is an n-th one of said 
mined phase reference. receivmg locations, and P„ is the complex detected 

49. The method of claim 48, further comprising the of the receiver located at and is expressed as: 
step of processing said detected complex values of the _ 

received sonic signals to determine phase velocity of 35 ^ u + # 

said multipole sonic waves in said borehole. -.^ ™ r 1 - • .i. 

50. nie method of claim 41, wherein step (d) com- . ^ The method of claim 53 further compnsmg the 
prises the step of calculating multipole sonic wave step of calculatmg the mulupole sonic wave renecUon 
phase velocity v and mulitpole sonic wave attenuation t*^""* °^ determined values for 
constant a from said detected complex values at two of ^ • ^ * , . ^ ^ 

said receiver locations and from the expressions: J*' »f ^bod of claim S3 further comprising the 

Step of calculating the multipole some wave phase ve- 
v(/)-Uiphi$e<p^2i,y))-ph«e(P^j2. locity V from the expression: 

45 v=(2if/)/Re»l Part (Jk) 

a= \n\mP^Z2, S^/P^zx, JWizx -22) 

56. The method of claim 53, further comprising the 
where the shear-wave detected-signal function is ex- step of calculating the multipole sonic wave attenuation 
pressed as: constant a from the expression: 

P^fi^^AifH-^t^^^^/^y^ ^ a=Im.giMryP«rt(k) 

and where P^is the detected-signal function; A is ampli- 57, The method of claim 41, wherein step (d) com- 

tude as a function of frequency and is independent of prises the steps of 

spacing between said emitting location and said receiver (a) defining the detected signal P(z) at a fixed fre- 
locations; z is the spacing between said emitting loca- quency to be 
don and a receiver location; f is the discrete frequency 
at which the multipole sonic wave energy is propa- 
gated; and Z] and zi are respective said receiver loca- /^j) » 2 ^„ e"^'^; 
tions. m=ai 

51. The method of claim 41. wherein step (d) com- ^ 

prises the step of calculating multipole sonic wave 0>) dc^'^B a linear predictor h(z) so as to satisfy the 

phase velocity v and multipole sonic wave attenuation relationship 

constant a from said detected complex values at three 

of said receiver locations spaced at intervals Az and W©*W=o 

from the expression: ^ , ^ . , . . 

where the notation. Qz), means convolution m z; 

*s=2ff(r/v)+/a=(iit G)/ii&3) (c) defining H(k) as the Fourier transform of h(z) in 

the wave-number domain, where said definitions of 



10/09/2004, EAST Version: 1.4.1 



27 



5.077,697 



receiver response P(2) and said linear predictor 
h(z) can be satisfied only when H(k) is rero at 

(d) truncating spatial variables P(z) and h(z) to finite 
intervals of z; ^ 

(e) solving the linear equations of the relationship 



Ar)A(r)-0 



10 



for h(z), assuming h(0)=: 1; 
(0 transfonaing the function b(z) to H(k); and 
(g) calculating values of k for the equation H(k)^0 to 

obtain estimated values for k^- 
58. The method of claim 57, further comprising the IS 
step of calculating the multipole sonic wave phase ve* 
locity and multipole sonic wave attenuation from said 
values for km from the relationship 



20 



59. The method of claim 57, further comprising the 
step of estimating said amplitude functions. by least 
square fitting. 25 

60. The method of claim 41, wherein step (d) com- 
prises the steps of 

(a) calculating the function L(k) as 



where 



GCA) = 



; and 
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and where R=[P][P]* is the covariance matrix 
firom the measured data, Pnis the complex receiver ^ 
response of the receiver located at Zn, N is the 
number of receivers in the array, and the super- 
script * denotes conjugate transpose of the matrix 
or the vector; and 

(b) determining values of km by searching for the 55 
locations of the local maxima of L(k). 

61. The method of claim 41, wherein step (d) com- 
prises the steps of 

(a) calculating the fimction p(k)=G»(k)R G(k), 
where 



60 
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and where R=[P][P]* is the covariance matrix 
from the measured data, P^ is the complex receiver 
response of the receiver located at z„, N is the 
number of receivers in the array, and the super- 
script * denotes conjugate transpose of the matrix 
or the vector; and 
(b) determining values of wavenumbers km by search- 
ing for the locations of the maTiTna of L(k). 

62. The method of claim 41, further comprising the 
step of deriving a radial shear velocity profile for a zone 
near said borehole. 

63. The method of claim 62, wherein said step of 
deriving a radial shear velocity profile comprises: 

(a) estimating a flexural wave velocity versus fre- 
quency profile in a zone surrounding said borehole 
from said detected complex values; 

(b) calculating a theoretical multipole wave velocity 
versus frequency and the shear velocity profile for 
said zone; and 

(c) searching for a flexural wave velocity versus fre- 
quency profile for said zone by curve-fitting until 
the diffcreocc between said estimated profile and 
said calculated profile is minimized. 

64. The method of claim 41, further comprising the 
step of determining shear velocity in a zone near said 
borehole by selecting the multipole sonic wave phase 
velocity at the lowest said discrete frequency for which 
a reasonable signal-to-noise ratio is obtained. 

65. The method of claim 41, further comprising the 
step of deriving a shear velocity value for a zone near 
said borehole. 

66. The method of claim 65, wherein said step of 
deriving a shear velocity value comprises: 

(a) estimating a multipole sonic wave velocity versus 
frequency profile in a zone surrounding said bore- 
hole from said detected complex values; 

(b) calculating a theoretical multipole sonic wave 
velocity versus frequency profile for said zone; and 

(c) searching for a shear wave value for said zone by 
curve-fitting until the difference between said esti- 
mated profile and said calculated profile is mini- 
mized. 

67. The method of claim 66, wherein said step of 
searching for a shear wave value comprises the step of 
finding the shear velocity Vcj^co, such that: 



/=l ^ C/3 - minimum, 

where vy? is the measured flexural- or screw-wave or 
other higher-order multipole wave velocity, v/f is the 
theoretically-calculated flexural- or screw-wave or 
other higher-order muldpole wave velocity, and >y/is a 
weighting function which can be chosen to be one or a 
function of wave amplitude at the i-th frequency. 

68. The method of claim 66, wherein said step of 
searching for a shear wave value comprises the step of 
finding the shear velocity vMea^such that: 
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niax \vn (fit - vn* is minimmn, 

. . . ^ 

where vjj is the measured flcxural- or screw-wave or 
other higher-order multipole wave velocity, vyy* is the 
theoretically-calculated flexural- or screw-wave or 
other higherorder multipole wave velocity, and w/is a 
weighting function which can be chosen to be one or a 
function of wave amplitude at the i-th frequency, 

69. Apparatus for borehole logging, comprising: 

(a) multipole source means for emitting sonic energy 
at a l^t location in a borehole, said sonic energy 
being of at least one discrete frequency which in- 
duces propagation of multipole sonic waves in the 
wall of said borehole; 

(b) multipole detector means for receiving sonic en- 
ergy of said multipole sonic waves at multiple re- 
ceiver locations in said borehole to produce a re- 
ceived sonic signal for each location, said multiple 
receiver locations being spaced apart from one 
another and from said first location; 

(c) means for detecting, for each said receiver loca- 
tion, the complex values of the received sonic sig- 
nal; and 

(d) means for processing said detected complex val- 
ues of the received sonic signals to determine atten- 
uation of said multipole sonic waves in said bore- 
hole as a function of frequency. 

70. The apparatus of claim 69, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
continuously emit sonic energy at said discrete fre- 
quency. 

71. The apparatus of claim 69, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
intermittently emit sonic energy at said discrete fre- 
quency. 

72. The apparatus of claim 69, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 
simultaneously emit sonic energy at a plurality of dis- 45 
Crete frequencies. 

73. The apparatus of claim 69, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer with 

a repetitive signal such that said source transducer emits 50 
sonic energy at a plurality of discrete frequencies below . 
6 KHz. 

74. The i^paratus of claim 69, wherein said emitting 
means comprises a multipole source transducer and 
means for driving said multipole source transducer to 55 
emit sonic energy at in sequential bursts at a variety of 
discrete frequencies, and wherein, for each said burst, 
said multipole detector means detects for each said 
receiver location amplitude of the received sonic en- 
ergy of said multipole sonic waves and phase of the 60 
received sonic energy of said multipole sonic waves 
relative to a phase reference. 

75. The apparatus of claim 69, wherein said multipole 
detector means detects for each said receiver location 
amplitude of the received sonic signal of said multipole 65 
sonic waves and phase of the received sonic signal of 
said multipole sonic waves relative to a phase reference 
derived from said emitted sonic energy. 



76. The apparatus of claim 69, further comprising 
means for processing said detected complex values of 
the received sonic signals to determine propagation 
velocity of said multipole sonic waves in said borehole. 

77. The apparatus of claim 69, further comprising 
means for processing said detected complex values of 
the received sonic signals to determine propagation 
velocity of said multipole sonic waves in said borehole 
as a function of frequency. 

78. The apparatus of claim 77, further comprising 
means for deriving a radial shear velocity profile for a 
zone near said borehole. 

79. The apparatus of claim 77, further comprising 
means for deriving a value of shear velocity for a zone 
near said borehole. 

80. A method of borehole logging, comprising the 
steps of: 

(a) einitting sonic energy from a multipole source at a 
first location in a borehole, said sonic energy being 
of at least one discrete frequency which induces 
propagation of multipole sonic waves in the wall of 
said borehole; 

(b) receiving sonic energy of said multipole sonic 
waves at multiple receiver locations in said bore- 
hole to produce a received sonic signal for each 
location, said multiple receiver locations being 
spaced apart from one another and from said first 
location; 

(c) for each said receiver location, detecting the com- 
plex values of the received sonic signal; and 

(d) calculating multipole sonic wave phase velocity v 
and multipole sonic wave attenuation constant a 
from said detected complex values at two of said 
receiver locations and from the expressions: 

/))]/2jrXzi-22) 
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where the shear-wave detected-signal function is 
expressed as: 

and where is the detected-signal function; A is 
amplitude as a function of frequency and is inde- 
pendent of spacing between said emitting location 
and said receiver locations; z is the spacing be- 
tween said emitting location and a receiver loca- 
tion; f is the discrete frequency at which the multi- 
pole sonic wave energy is propagated; and z\ and 
Z2 are respective said receiver locations. 
81; The method of claim 80. wherein step (a) com- 
prises continuously emitting said sonic energy at said 
discrete frequency. 

82. The method of claim 80, wherein step (a) com- 
prises intermittently emitting said sonic energy at said 
discrete frequency. 

83. The method of claim 80, wherein step (a) com- 
prises simultaneously emitting said sonic energy at a 
plurality of discrete frequencies. 

84. llie method of claim 80, wherein step (a) com- 
prises repetitively emitting said sonic energy at a plural- 
ity of discrete frequencies below 6 kHz. 

85. The method of claim 80, further comprising the 
step of repeating steps (a)-(c) for each discrete fre- 
quency. 
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96, The method of claim 80» further comprising the 
top of deriving a predetermined phase reference from 
said emitted sonic energy, and wherein said step of 
detecting the complex values of the received sonic sig- ^ 
na] comprises the step of detecting said complex values 
of the received sonic signal in relation to said predeter- 
mined phase reference. 

87. A method of borehole logging, comprising the 
steps of: 10 

(a) emitting sonic energy from a muldpole source at a 
first location in a borehole, said sonic energy being 
of at least one discrete frequency which induces 
propagation of muldpole sonic waves in the wall of 
said borehole; 

(b) receiving sonic energy of said multipole sonic 
waves at multiple receiver locations in said bore- 
hole to produce a. received sonic signal for each 
location, said multiple receiver locations being 20 
spaced apart from one another and from said first 
location; 



32 



(c) for each said receiver location, detecting the com- 
plex values of the received sonic signal; and 

(d) calculating multipole sonic wave phase velocity v 
and multipole sonic wave attenuation constant a 
from said detected complex values at three of said 
receiver locations spaced at intervals Az and from 
the expression: 

A=2irC//v)+id=(ln G)/(/Ar) 
where 



2P2 



and where k is the complex wave-nuMber of the 
received multipole sonic wave, f is the frequency of 
the received multipole sonic wave, and Pn is the 
complex detected signal of the receiver located at 
an n-th one of said receiving locations. 
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